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Abstract 

The recent observation of and B^ — > decay modes allows us 

to make a fresh isospin analysis of B ^ vrvr transitions. We find that current 
experimental data can impose some model-independent constraints on the 
parameter space of direct CP violation in B^ — > tt^tt" and Ba — > 
Furthermore, we establish a direct relationship between the weak phase a 
and the charge-averaged branching fractions and CP-violating asymmetries 
oi B ^ vrvr decays. 
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The major goal of B-mcson factories is to test the Kobayashi-Maskawa mechanism of 
CP violation within the standard model and to detect possible new sources of CP violation 
beyond the standard model. So far the CP- violating asymmetry in vs — > J/ipK^ 
decays has unambiguously been measured at KEK and SLAC [1] , and the experimental result 
is compatible very well with the standard-model expectation. Some preliminary evidence 
for CP violation in a number of other B decays has also been reported by BaBar and 
Belle Collaborations [2]. Further experiments will provide much more precise data on CP 
violation in the 5-meson system, from which one may cross-check the consistency of the 
Kobayashi-Maskawa picture and probe possible new physics. 

The charmless two-body nonleptonic decays B^ tt+tt", B^ — > n^n^ and B~ — > 7r°7r~, 
which can be related to one another via the isospin triangle [3], have been of great in- 
terest in B physics for a stringent test of the factorization hypothesis, a quantitative 
analysis of final-state interactions, and a clean determination of the CP-violating phase 
a = arg[—{VtdVfl) / {VudV*fj)] (an inner angle of the well-known unitarity triangle of the 
Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix [4]). Among these three decay 
modes, B^ h^tc^ is color-suppressed and should have the smallest branching fraction. 
The observation of this rare decay has recently been reported by BaBar [5] and Belle [6] 
Collaborations ^: 



'00 



_ ( (2.1 ± 0.6 ± 0.3) X 10-6 (BaBar) , , 

~ I (1.7 ± 0.6 ± 0.3) X 10-6 (Belle) , ^ ^ 

where the first error is statistical and the second one is systematic. The experimental value of 
Boo is a bit larger than the theoretical prediction in various QCD models [7]. In comparison, 
the observed branching fractions of B^ tt+tt- and B~ tt^tt- decays [8] 



g(^g^7r+7r-)+g(^°^7r+7r-) 
2 

' (4.7 ±0.6 ±0.2) X 10-6 (BaBar) 

< (4.4 ± 0.6 ± 0.3) X 10-6 (Belle) (2) 
(4.5 ± 1.4 ±0.5) X 10-6 (CLEO) 



and 

Bo± = 



_ Bj Bj ^ TT^TT^) + B{B- ^ TT^TT-) 

2 

' (5.5 ± 1.0 ±0.6) X 10-6 (BaBar) 

= < (5.3 ± 1.3 ± 0.5) X 10-6 (Belle) (3) 

(4.6 ± 1.8 ±0.7) X 10-6 (CLEO) 

are essentially consistent with the theoretical results [7]. Because Bqq and Bo± are all 
charge- averaged, they cannot directly be related to one another through the isospin relation. 



Note that the experimentally-reported branching fractions of 5 ^ tttt decays are all charge- 
averaged. 
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The purpose of this note is to make a fresh isospin analysis of three B tttt decays by 
using new experimental data and taking account of CP violation. We find that current data 
can impose some useful and model-independent constraints on the parameter space of direct 
CP violation in B^ — > tt+tt" and B^ — > 7r°7r° transitions. It is also possible to establish a 
direct relationship between the weak phase a and the charge- averaged branching fractions 
and CP-violating asymmetries of B ^ tttt decays. 

Let us define the direct CP- violating asymmetries between B^ 7r"'"7r~, B^ 
B~ — > 7r°7r~ and their CP-conjugate decays: 
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(4) 



B{B- ^ 7r%-) + B{Bt 
Taking account of Eqs. (l)-(4) , one can easily obtain the magnitude of each decay amphtude: 
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(6) 



where Tq and Tj_ denote the lifetimes of neutral and charged B mesons, respectively. The 
present experimental data yield k = Tq/t^ ~ 0.92 [4]. In Eqs. (5) and (6), we have neglected 
the tiny phase space difference between tt^tt^ and 7r"'"7r~ (or tt^tt^) states. It is obvious that 
the relevant decay amplitudes cannot be determined, unless the CP-violating asymmetries 
Ao and A± are all measured. 
Under isospin symmetry and in the neglect of electroweak penguin contributions [9], the 



amplitudes of 5° — >• tt"*"?! , B\ 
the complex plane: 



7r°7r° and B^ tt^tt decays form an isospin triangle in 



(7) 
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Similarly, 



(8) 



In this approximation, the magnitudes of A{Bl^ tt^tt") and A{B^ tt^tt^) are identical 
to each other [3]. Thus the CP- violating asymmetry ^o± vanishes. What we are concerned 
about is the relative phase between the amplitudes of — > n'^n^ (or B^ — > tt+tt") and 
7r°7r° (or 5° ^ n\°) decays: 



= arg 



ip = arg 



MB', 



MB', 



71^71 



A(50 ^ vrOvrO) 



TT+TT ) 



(9) 



At the tree level 
effects in B, 



(f = (f) holds. A difference between and ip measures the QCD penguin 
TT^TT" and Bd tt^tt^ transitions. With the help of Eqs. (5)-(8), we find 



COS0 = 



cos (fi — 



2kB, 



'0± 



2i3oo (1 + Aoo) 



2^2B+-Boo (1 + A+-) (1 + Ao) 
2kBo± - g+_ (1 - A-) - 2Boo (1 - Ao) 
2^2B+_Boo (1 - A-) (1 - Ao) 



(10) 



where A± = has been taken into account. One can sec that Lp = (f) would hold, if 
both A+- and Ao were vanishing. Hence the difference between and (f results from the 
penguin-induced CP violation in B, —>■ tt+tt" and B, 7r°7r° decays. 

To get a ball-park feeling of the discrepancy between and (/?, we make use of the 
QCD factorization to estimate the ratio of A{B'^ tt+tt") to A{B^ — >■ tt^tt^) and that of 



MB', 



TV ' TT 



to A(i?o 



TT ' TT 



[10]. The results are 
(ai + a2 + ag^^) RbC^ + {al + Og^^) 



(02 



I - <e.) RbC-'^ - (at + ageo 

(ai + < + <e^) i?fee+^^ + {at + age^r) 

A(fiO -> ttOttO) (as - a| - <G) Rbe+'^ - {ai + a^n) 



MB', 



71^71 



(11) 



where Rb = \Vu,V*f,\/\Vc,V*f,\ and 7 = £Lrg[-{Vu,V*f,)/{Vc,VX)] are two parameters of the 
CKM unitarity triangle [11]; and a^'^ denote the effective Wilson coefficients [12]; 

is associated with the penguin operators (^5,6 [13] and can be given as 



2mt 



(nib - my,) {niu + m,) 



under isospin symmetry. For illustration, we typically adopt Rb 
as well as ai ^ 1.038 + O.OlSi, 02 ~ 0.082 - 0.080i, < fa - 



(12) 

0.37 and 7 « 65° [11] 
0.029 - 0.015i and al ^ 



-0.034 — 0.008i at the scale n = rrib [12]. The formally power-suppressed QCD coefficients 
q''^ can be neglected in the heavy quark limit [12], because of ~ Then we arrive at 



MB', 



7r+7r-)/A(E0 ^ ttOttO) ^ 7Ae-'^^ and ^1(5° ^ 7r+7r-)/A(E0 



4.3e+'^4°. 
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As a result, (f) ~ —17° and ip ^ 44° are obtained. The big difference 99 — ^ 61° implies 
the existence of significant QCD penguin effects and large direct CP violation in B^^ tttt 
decays. It should be noted that and ip are neither pure weak phases nor pure strong phases, 
since the treel and penguin amplitudes of each —>■ tttt decay mode involve different weak 
and strong phases. It should also be noted that the numerical results for and (f are just 
illustrative, because possible final-state rescattering effects have not been taken into account 
in Eq. (11). 

Once the CP- violating asymmetries and ^00 ^^I'c measured at B factories, it will be 
possible to determine (j) and (p from Eq. (10) independently of the specific QCD models. A 
hint of A+- 7^ (of ~ 2a significance) has recently been reported by the Belle Collaboration, 
but the BaBar Collaboration's preliminary result for A+- is consistent with zero [2]. There is 
no experimental information about ^00 at present. In this situation, we show the numerical 
correlation between 4> and p in Fig. 1(a) and that between A+- and ^00 in Fig- 1(b), where 
current data fisted in Eqs. (l)-(3) have been used. The points in Fig. 1(a) are generated 
from scanning the error bars of Bqq and Bo±. The fight region in Fig. 1(b) is generated 
in a similar way and corresponds to the allowed ranges of cos0 and cosp in Fig. 1(a), while 
the dark region in Fig. 1(b) results from inputting the central values of B+-, Bqo and Bo±- 
Some comments and discussions are in order. 

(1) Current experimental results of B^_, Bqq and Bo± can impose some limited but 
model-independent constraints on the values of (0, </?) and (.4+_,.4oo)- For example, the 
possibilities cos0 = cosp = and = |.Aoo| = 1 are not allowed. As a consequence of 
|^+_| < 1 and I ^00 1 < 1, the regions for cos0 < —0.2 and cos(/7 < —02 have completely 
been excluded by the present experimental data. 

(2) The error bars of Bqq and Bo± influence the parameter space of (0,93) or 
(^+_,^oo) to a very limited extent. If A+- is measured, one can get a roughly allowed 
range of .4oo from Fig. 1(b). For instance, A+- > 0.6 and .4oo > 0.6 cannot simultaneously 
hold. Nor can A^ < —0.6 and ^00 ^ —0.6 simultaneously hold. 

(3) As the tree-level amplitude of n^ii^ or vr^vr^ is color-suppressed, its 
magnitude might be comparable with the size of the corresponding penguin amplitude. 
Then |^oo| is expected to be larger than |^+_| in the QCD factorization. However, the 
elastic rescattering effects in the final states of S — > tttt transitions could easily spoil this 
naive expectation. It is therefore useful to analyze the relevant experimental data in a 
model-independent way. 

(4) The validity of our numerical results relies on the validity of the isospin relations in 
Eqs. (7) and (8), which hold in the assumption of negligible electroweak penguin effects. 
The electroweak penguin contributions to S — > tttt transitions are generally expected to 
be insignificant [9]. This expectation would be problematic or incorrect, if Ao± 7^ were 
experimentally established. It is also worth mentioning that final-state interactions in 5 — >■ 
Tin decays include both elastic nn ^ tttt rescattering and some possible inelastic rescattering 
effects ^. Whether inelastic final-state interactions are negfigibly small or not remains an 



^The amplitudes of tt+tt , 5° vr'^Tr^ and tt^tt decay modes (or their CP-conjugate 

processes) may still form an isospin triangle in the complex plane, even if the inelastic tttt ^ DD 
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open question [16]. To answer this question requires more precise data on both the branching 
fractions and the CP- violating asymmetries oi B ~> im decays. 

We proceed to discuss the extraction of the CP-violating phase a with the help of the 
isospin triangles in Eqs. (7) and (8). As shown in Ref. [3], a depends on two CP- violating 
observables which arise from the interplay of decay and B^-B^ mixing in B, — > tt'^tt" and 
Bd — > 7r°7r° decays: 



X+- = 



Xoo 



Im 


Q 


MBl - 




P 


MB'a - 


-> 7r+7r~ 


Im 
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|i?|sin[2(a + e)] , 
|i?|sin[2(a + e)] , 



(13) 



where 



R = \R\e 



2ie 



1-r ' 
2 + f 



2 + r ' 



(14) 



with r = |r|e'^ = A0/A2 and f = |f|e*^ = Ao/ A2 being the ratios of / = and 1 — 2 isospin 
amplitudes in B^ — > tt'^'tt" (or 7r°7r°) and B'^ n'^n' (or 7r°7r°) decays. Following Ref. [17], 
we obtain 



and 



where 



9 — ± arccos 



3(a + 6) - 2 , 

66 — 3a 



4J3(a + 6) -2 



9 — ± arccos 



3(a + 6) -2 , 
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(15) 



(16) 



(17) 



rescattering effects are taken into account [14]. In this specific case, the analytical result obtained 
in Eq. (10) remains valid. It is even expected that inelastic final-state interactions could account 
for the existing discrepancy between experimental data and QCD models for D ^ tttt and B ^ nir 
decays [15]. 
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a= ^^^.^^ ^ , 

(l + Ao). (18) 

Clearly r (or f) can be determined up to a twofold ambiguity in the sign of its phase, if both 
A+- and Aoo are measured. One can in turn determine R (or R), both its magnitude and 
its phase. It is then possible to extract the weak phase a from Eq. (13) through the time- 
dependent measurement of x+- and (or) xoo in Bd — > tt'^tt' and (or) Bd 7r°7r° decays. On 
the T{4S) resonance, the time-dependent rate asymmetry between B^{t) and B^{t) decays 
into a CP eigenstate / reads as 



A/(t) 



mit)^f]+rmt)^f] 

\Pf\' - 1 




cos(AMd ■ t) + , , / sin(AMrf ■ t) , (19) 



where pj = A{B^ f)/A{B^ — > /) and tiny CP violation in B^-B^ mixing has been 
neglected (i.e., \q/p\ ~ 1) [18]. For / = tt+tt" and 7r°7r°, we take account of Eqs. (4)-(6) 
and (13) and obtain 

A,+,- (t) = A+- cos(AMd ■t) + {l- A+-)x+- sm{AMd ■ t) , 

A^o^o (t) = Aoo cos(AMd -t) + (1 - Ao)xoo sin( AM^ • t) . (20) 

Clearly both A+^ (or Aoo) and x+- (oi' Xoo) can be determined from measuring the time 
distribution of A^r+Tr- (or A^o^o) at asymmetric B factories. Compared with Refs. [3,17], the 
formulas obtained in Eqs. (13)-(18) and (20) establish a more direct relationship between 
a and the observables Bqq, Bo±, A+- and Aoo- Thus the usefulness of our results is 
noteworthy. 

To summarize, we have presented a fresh isospin analysis of rare B tttt decays by 
taking account of the fact that the experimentally-reported branching fractions are charge- 
averaged and large CP violation may exist in them. We find that current experimental data 
can impose some useful and model-independent constraints on the parameter space of direct 
CP violation in Bd tt+tt" and Bd — > 7r°7r° transitions. A straightforward relationship 
between the CP-violating phase a and the measurables of B ^ tttt decays has also been 
established. More precise measurements of such charmless B decays will allow us to test the 
consistency of the Kobayashi-Maskawa mechanism of CP violation and probe possible new 
physics beyond the standard model. 

The author likes to thank Y.F. Zhou for useful discussions and partial involvement at the 
very early stage of this work (in 2002). This research was supported in part by the National 
Nature Science Foundation of China. 



7 



REFERENCES 



[1] BaBar Collaboration, B. Aubert et al, Phys. Rev. Lett. 87, 091801 (2001); Belle Col- 
laboration, K. Abe et al, Phys. Rev. Lett. 87, 091802 (2001). 

[2] H. Yamamoto, plenary talk given at the International Europhysics Conference on High 
Energy Physics, July 2003, Aachen, Germany; T. Browder, talk given at the 21st In- 
ternational Symposium on Lepton and Photon Interactions at High Energies (LP 03), 
August 2003, Batavia, Illinois, USA; H. Jawahery, talk given at the same conference 
(LP 03). 

[3] M. Gronau and D. London, Phys. Rev. Lett. 65, 3381 (1990); Y. Grossman and H.R. 

Quinn, Phys. Rev. D 58, 017504 (1998). 
[4] Particle Data Group, K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 
[5] BaBar Collaboration, B. Aubert et a/., hep-ex/0308012. 
[6] Belle Collaboration, K. Abe et al, hep-ex/0308040. 

[7] For recent reviews with extensive references, see: H.N. Li, hep-ph/0303116; M. Beneke 

and M. Neubert, hep-ph/0308039. 
[8] J. Fry, talk given at the 21st International Symposium on Lepton and Photon Interac- 
tions at High Energies (LP 03), August 2003, Batavia, Illinois, USA. 
[9] N.G. Deshpande and X.G. He, Phys. Rev. Lett. 74, 26 (1995); M. Gronau, O.F. 
Hernandez, D. London, and J.L. Rosner, Phys. Rev. D 52, 6374 (1995); R. Fleischer, 
Phys. Rept. 370, 537 (2002). 
[10] Z.Z. Xing, Phys. Lett. B 493, 301 (2000). 
[11] A.J. Buras, hep-ph/0307203. 

[12] M. Beneke, G. Buchalla, M. Neubert, and C.T. Sachrajda, Phys. Rev. Lett. 83, 1914 

(1999); Nucl. Phys. B 591, 313 (2000). 
[13] G. Buchalla, A.J. Buras, and M.E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 
[14] C. Hamzaoui and Z.Z. Xing, Phys. Lett. B 360, 131 (1995). 
[15] H.Y. Cheng, Eur. Phys. J. C 26, 551 (2003). 

[16] J.F. Donoghue, E. Golowich, A.A. Petrov, and J.M. Soares, Phys. Rev. Lett. 77, 2178 

(1996); Z.Z. Xing, Phys. Rev. D 53, 2847 (1996). 
[17] D. Du and Z.Z. Xing, Phys. Rev. D 47, 2825 (1993). 
[18] See, e.g., Z.Z. Xing, Phys. Rev. D 53, 204 (1996). 



8 



FIGURES 



cos (p 




COS 



0.5 




-0.5 - 



00 



FIG. 1. Constraints of current experimental data on (a) the parameter space of cos ^ and cos 
and (b) the parameter space of ^oo and A+-. Note that the dark region of (b) corresponds to the 
central values of i3oo and Bq± given in Eqs. (l)-(3). 
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